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Abstract involvement has been at the forefront of these 

developments. 


A review of analytical methods and associated 
tools for active controls analysis and design 
problems is presented. Approaches employed to 
develop mathematical models suitable for control 
system analysis and/or design are discussed. 
Significant efforts have been expended to develop 
tools to. generate the models from the standpoint 
of control system designers' needs and develop the 
tools necessary to analyze and design active 
control systems. Representative examples of these 
tools are discussed. Examples where results from 
the methods and tools have been compared with 
experimental data are also presented. Finally, a 
perspective on future trends in analysis and 
design methods is presented. 


Introduction 


Active controls technology offers the 
potential for realizing economic and performance 
benefits from aircraft configurations that would 
be unacceptable by traditional design standards. 
Through interdiscipl inary efforts, especially 
among controls specialists and aeroelasticians, a 
new concept is evolving for aircraft design in 
which both aircraft stability and structural 
integrity are dependent on the operation of an 
active control system. Because active controls 
analyses span the traditional disciplines of 
structural dynamics, aerodynamics (both steady and 
unsteady), propulsion (propulsion will not be 
considered in this paper), and control theory, a 
common format for the math models is required. 1 
The integration of these traditional disciplines 
is maturing into a new discipline termed by many 
as aeroservoelasticity. 

The analytical methods (both analysis and 
design) have been evolving for some time. 2-9 in 
design, the problem is the synthesis of a control 
law that regulates the dynamics of the system such 
that measures of performance (e.g. , stability, 
loads, transient response, etc.) are acceptable. 

In analysis, the control law already exists and 
analytical methods are employed to assess the 
performance of the overall controlled system. 

To take full advantage of active control 
technology, control law synthesis and analysis 
must be an integral part of the aircraft design 
process. This requires efficient methods and 
accompanying tools that will enable the active 
controls designer routinely to synthesize and 
analyze complex control systems. A significant 
amount of work has been performed to develop 
analysis and synthesis methods and tools for this 
task. Much of the NASA active controls research 


Since the validation of analytical methods 
is, as always, very important, several experi- 
mental studies have been performed to date. 10 - 1 ? 
These experimental studies have ranged from wind- 
tunnel tests to subscale flight tests to full- 
scale flight tests. The majority of the 
wind-tunnel tests have emphasized flutter 
suppression whereas the majority of the flight 
tests have emphasized other active control 
concepts such as relaxed static stability, 
maneuver and gust load alleviation, etc. 

The purpose of this paper is to describe the 
development and validation of analytical methods 
for application to active controls technology. 
Although various analytical methods will be 
mentioned, emphasis is given to the methods 
developed at the NASA Langley Research Center. A 
few representative examples of studies to 
correlate analysis and experimental data are 
presented. Finally, some future trends in 
analytical methods will be discussed. 


Nomenclature 


Aj real aerodynamic approximation 

matrices 

A » B u» B w q P^nt matrices 
A c ,B c ,CH controller matrices 
A a augmented plant matrix 

B a augmented input matrix 

b reference length 

C output matrix 

C co control law optimum gain matrix 

D v viscous damping matrix 

E expected value 

F(s) response transfer function matrix 
H(s) modal transfer function matrix 

i fA 

J cost function 

K$ generalized stiffness matrix 

k reduced frequence, bw/V 

L Kalman filter gain matrix 

M generalized mass matrix 

N 0 average no. of positive zero 

crossing 

p nondimensional Laplace variable, 

bs/V 

Q force vector function 

Ql,Q 2 weighting matrices 

q dynamic pressure 

noise intensity matrix 
s Laplace operator 

t t i me 

u plant input vector 

U system input function vector 
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velocity 
gust vector 

measurement noise vector 
vector of state variables 
augmented state covariance matrix 
vector of controller states 
output vector 

vertical displacement at point 

(x,y) at time, t 

aerodynamic lag terms 

control surface deflection 

actuator output 

actuator input 

jth mode shape 

gust input power spectrum 

phase angle 

noise vector 

rms value 

circular frequency 

jth generalized coordinate vector 

damping coefficient 


Equations of Motion 


Modeling the dynamics of an actively 
controlled flexible vehicle is a multidisciplinary 
task. Structural dynamics, unsteady aerodynamics, 
and control elements are required. Structural 
dynamics, unsteady aerodynamics, and control 
system interactions can be represented by 
second-order frequency domain or first-order 
state-space (time domain) equations. A flow chart 
which illustrates the development of these 
equations is shown in figure 1. A description of 
the development of the mathematical models 
follows. 


Structural Dynamics 

The" structural" dynamics for a flexible 
airplane can be formulated using Lagrange’s 
equations. 18 It is normally assumed that the 
dynamics can be represented by a finite number of 
modes. By the method of separation of variables, 
the motion is assumed to be the product of a shape 
function and a time function. 

N 

*(x,y,t) =!E ♦j(x,y)£j{t) (l) 

3=1 


If all of the structural damping in the aircraft 
is assumed to be viscous in nature, then the 
equations of motion can be written as 


for stable systems, for computing a statistical 
measure of the response to random inputs. 

However, Q(p) is needed for investigations of 
system stability. This point will be discussed 
further in a later section entitled “Stability." 

The concept of analytic continuation 6 can 
be employed to develop an approximation to the 
aerodynamic forces for arbitrary motion (i.e., a 
function of p) in terms of known aerodynamic 
forces for oscillatory data. 1 * 4 * 6 * 21 The most 
widely used form is 

L A i+2 P 

A 0 + A x p + A 2 p2 + l _____ (3) 

3=1 Bj+P 

The real coefficients are determined, subject to 
any imposed constraints (e.g., Aq associated 
with rigid-body linear displacement should be 
zero), such that the error between the approxi- 
mation and the known tabular data is minimized 
(e.g., in a least squares sense) 1 * 4 * 21 

The motivation for the rational approximation 
of equation 3 is to enable the equations of motion 
to be transformed to a set of constant-coefficent 
first-order differential equations that are 
amenable to analysis by efficient linear systems 
techniques. Such a transformation will be defined 
in a section entitled “State Space Equations of 
Motion." 


Control System 

” Control system equations are normally 
described by transfer functions. Standard 
transformation techniques can be employed to 
transform the transfer function representation to 
a state-space representation. The general 
relationship between actuator outputs and inputs 
can be expressed as 

*A = {T a (s)K(s)T s (s)C(s) + T he (s)H e } 5 + T A (s)6 com 

( 4 ) 

The definition of these matrices will be given in 
the next section. 


Overall System 

Bycombini ng the structural dynamics, 
unsteady aerodynamics, and control system 
representations, the equations of motion In the 
frequency domain (see figure 2a for a block 
diagram) can be written as 


Ml + 0 V | - K s £ = Q 


(2) {Ms 2 + D y s+K s +qQ(^)' - F D d HS (s )C (s ) 


The forces on the right hand side of equation 2 

are a result of aerodynamic forces due to aircraft + T HE (s ) H£>} C = w- Q w + F n T A^ s ^com 

motion, turbulence, and control surface motions * 99 m (5) 

and forces due to actuator motions. 

where 


Unsteady Aerodynamic Forces D v 

Unsteady aerodynami c forces are normally q 

computed as tabular functions of Mach number and V 

the complex variable p=*sb/V. Programs currently b 

available for production-generation of aerodynamic p=sb/V 

forces can only perform the computations for p* 

/-l b^/V. 19 * 20 This is precisely the form 
required for frequency response computations or. 


viscous damping coefficient 
dynamic pressure 
aircraft speed 
reference length 
nondtmensional Laplace variable 



Q generalized aerodynamic force 

C(s) input to sensor i resulting from 

motion in the jth generalized 
coordinate 

Hg; N^xNp matrix of modal 

coefficients relating actuator 
displacement to generalized 
coordinates 

T HS { s ) f^xN s matrix relating actuator 

hinge moment outputs to sensor 
inputs. This matrix contains 
sensor, control logic, and actuator 
dynamics. 

Tre( s ) matrix relating actuator 

hinge moment outputs to actuator 
displacements. 

Tft(s) diagonal matrix of 

transfer functions relating actuator 
hinge moment outputs to actuator 
displacements 

T$(s) N s xN s diagonal matrix of 

transfer functions relating sensor 
outputs to sensor inputs 
Fp f^xNx ntetn* of modal 

coefficients converting hinge moment 
outputs to generalized forces 
6 com commanded (e.g., pilot) inputs to 

actuators 

a ^ outputs from actuators 

Wg gust velocity input 

An alternate form of the equations often used 
for simplicity neglects aerodynamic hinge moments 
and hinge moments due to inertial coupling between 
control surface and basic wing degrees of freedom 
and also assumes infinite stiffness of the backup 
structure. In this case control surface 
deflections are given by 

6 - T$(s)C{$)£ + fA( s )^com (6) 

The dimension of 6 is N 5 xl and T 6 Is a 
matrix of transfer functions relating control 
deflections to sensor inputs. The resulting 
equations can be expressed as 



(7) 


State Space Equations of Motion . Under the 
a ssumpt i on that the approximation TrT equat i on (3) 
is valid for points off the imaginary axis, linear 
time-invariant equations of motion can be written 
in the form (see figure 2b for a block diagram) 


X = AX + B U U + ^WgWg 


y = CX + Wjn 


( 8 ) 


where 
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u com 
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n^ x n^ states associated 

with the unsteady aerodynamic 
force representation 
states representing sensor 
dynamics transfer functions 
states representing actuator 
transfer functions 
gust states representing a 
filter which converts white 
noise input to a gust velocity 
having a power spectral density 
approximating atmospheric 
turbulence 

commanded control inputs 

uncorrelated zero mean, 

white noise processes with 

intensities Rg and ^ 

respectively 

states representing the 

control logic 

matrices in state space 

representation of the control 

logic 


Analysis Methods 


Analysis consists of solving the equations of 
motion for quantities that provide information on 
the performance of the active control system. An 
actively-controlled flexible aircraft may require 
analysis of many different measures of 
performance. Figure 3 shows plots typical of 
those used to assess the performance of an active 
control system. All of these plots will probably 
be examined during the design of an active control 
system. 


The analysis methods employed for an 
actively-controlled flexible aircraft can be 
divided into two major categories. The first 
category involves classical methods defined in 
this paper where the frequency -domain formulation 
of the equations of motion is used. The second 
category involves modern methods defined in this 
paper where the state-space formulation of the 
equations of motion is used. All of the 
quantities that are needed to assess the 
performance of an active control system can be 
found using either classical or modern methods. 
However, one method may be better suited to 
evaluate certain quantities. Techniques used in 
both categories of analysis methods will be 
described. 

Classical Analysis Methods 

As stated previously , classical analysis 
methods normally utilize the frequency domain 
formulation of the equations of motion. 

Equation 5 can be written as 


X c = + B C Y 

u = C c x c + u com 


(9) 


H(s)£ - F(s)U 


( 10 ) 
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The O represents a closed-loop quantity* 


the singular value computations is described in 
reference 25. 


Stability. The unsteady aerodynamic force 
matrices are currently only calculated at s^iw. 
(Reference 22 describes an approach which would 
allow direct computation of Q(s).) The availabil- 
ity of only Q(ioj) creates a problem in solving for 
the eigenvalues of H(s). Two approximations are 
commonly employed to circumvent this 
difficulty* 4 * 2 ^ 

One approximation, valid for lowly damped 
eigenvalues, is to assume 


Stochastic response . Classical analyses are 
typically per formed to examine a statistical 
characterization of loads due to turbulence at 
critical locations in an aeroelastic struc- 
ture. 26 * 27 Historical data have been used to 
establish certain guidelines that, if followed, 
result in a safe structure. With the advent of 
active controls, these guidelines have been 
applied for closed-loop cases. The guidelines 
have been extended, somewhat arbitrarily, to 
establish statistical control power criteria. 2 ^ 


Q(s) - Q(0+iu>) * Q(i Vk/b ) (11) 

and approximate eigenvalues are found by using the 
p-k method. 23 When the p-k method is employed 
to solve 


jTf(s)| = det Tf(s) ■« 0 

Q(ik) is a nonlinear, tabular function of k. 
Consequently , solutions can only be obtained by an 
iterative process. 

The other approach is to make a rational 
s -plane approximation 

ft(s) - Q(s) 

as in equation 3. Eigenvalue determination where 
one has Q(s)=Q(Vp/b) is commonly termed the p-p 
method. When Q(s) is employed, one may transform 
the equations of motion as in equation 8 into a 
standard eigenvalue form. Then the iterative p-p 
procedure and standard eigenvalue (state-space) 
methods produce identical results. 


If the gust input is a Gaussian random 
process, then two quantities, a (the rms value) 
and N 0 (the average number of positive zero 
crossings), suffice to describe the output 
response to the turbulence excitation. These 
quantitie: 

a- > 

2 "V 

"» ' 55T 

W J 

where fs the gust input power spectrum. 

Deterministic response . As stated before, 
the unsteady aerodynamic matrices are calculated 
only for s=iu> (steady-state oscillatory motion) 
and the equations of motion are written in the 
frequency domain. Fourier transform techniques 
provide a method by which deterministic response 
cart be obtained from the solution of the frequency 
domain formulation of the equations of motion. 


can be written as 

2 T 1/2 


'/ 

t . 0 

It 


% ( la) ) 
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$ w (a>)da> > 

9 J 
2 1 1/2 
$ w ^(u))du) \ 


(13) 


Figure 4 (from ref. 24) shows elastic mode 
eigenvalue loci as a function of feedback gain for 
a single-input-single-output control law. Curves 
are shown that were obtained by the p-k and the 
p-p methods. Differences between the loci found 
by the two methods increase as damping ratio 
increases* The slight difference at the flutter 
point is due to the fact that Q(lu) did not fit 
the Q(i<a) data precisely at the flutter 
frequency. 


Deterministic response analysis by Fourier 
transform techniques can be separated into three 
phases. First, the forcing function (defined in 
the time domain) is transformed into the frequency 
domain. Second, the responses are computed in the 
frequency domain. Third, the responses (in the 
frequency domain) are transformed back to the time 
domain. Reference 29 gives a complete description 
of this method to calculate deterministic 
responses. 


Robustness. The plant and controller 
i n pu t / ou tp u ifTFa n s f e r matrices can be written as 


V - T(iai)C(iai)H 1 (ioj)F(iaj)U - Gg(iaj)6 A +G w w g 
6 a = K(iu>)Y + T A (iu))6 com 


( 12 ) 


Classical analysis techniques have, in the 
past, relied heavily upon Bode plots of particular 
elements of % and K and also upon Hyquist plots 
of elements of the loop transfer matrix %K. 


The classical formulation is also particular- 
ly well suited for analysis of the eigenvalues and 
singular values of the loop transfer and return 
difference matrices. In the classical formulation 
no s-pl ane approximation (with accompanying addi- 
tional states) is needed and experimental frequen- 
cy response data for dynamic elements such as 
actuators and sensors can be utilized directly 
without adding additional states to represent 
them. Use of the classical approach to perform 


Modern Methods 

The commonly conducted analyses using the 
state-space representation include block- 
diagonal izati on , controllability and observ- 
ability, eigenvalue/eigenvector, covariance, fre- 
quency response, time response, and robustness. 

The matrix A (eq„ 8) is normally a highly 
coupled matrix. The coupling is due to presence 
of the aerodynamics and actuator dynamics. A 
linear transformation can be used9 to transform 
the system to block diagonal form. The diagonal 
blocks are composed of the real and imaginary 
parts of the eigenvalues of each mode. The 
transformed B and C matrices provide information 
regarding the degree of controllability and 
observabi 1 ity of each state. 

The characteristic roots of the system can be 
obtained by computing the eigenvalues of the 
matrix A and may be plotted for varying dynamic 
pressures or other parameters to provide root loci 
or parametric stability information., No 
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determinant Iteration Is required (unlike the p-k 
or p-p methods). 

For a stable system the root mean square 
(rms) response to a stationary white noise Input 
n(t) Is obtained using covariance analysls.30 
The computation Involves solution of a steady- 
state Lyapunov equation 


Va 


+ X a A a 


+ B R B R T = 
ana 


(14) 


where X a is the covariance matrix of the 
augmented state vector } and is the noise 

intensity matrix. The square root of the diagonal 
terms of the X a matrix are the rms values of 
each state. In classical analysis, the rms values 
are obtained by numerically integrating the 
power-spectral density over a finite-frequency 
domain. 


The time response to a deterministic input 
u(t) is given by 

X(t) = e A ^ t_t °^X(t 0 ) + f e A ^ t * T ^Bu(T)dT 
To 

The computation of the state transition matrix 
e^ t can be facilitated by block diagonal trans- 
formation of the matrix A. For certain u(t) time 
histories, the time Integration can be performed 
analytically. 


decoupled. 34 stable nominal closed -loop systems 
can be obtained using these approaches. However, 
additional computations are required to quantify 
robustness characteristics. Tests have been 
Identified 3 * which Involve both principal gains 
(singular values) and principal phases of the loop 
transfer matrix. These tests are less conserv- 
ative than robustness tests based solely upon 
singular values. 

Recently, classical (input/output) approaches 
which utilize fractional representation theory 36 
have been developed which allow a parametric defi- 
nition of the class of all stabilizing controllers 
for a given plant. Thus one can be assured while 
varying the free controller parameters to seek a 
more desirable subset, that the stability of the 
nominal closed- loop 'system is preserved. 

Modern methods . The most popular modern con- 
trol theory method is the Linear Quadratic 
Gaussian (LQG) method. 37 Because of the high 
order of the model for a flexible aircraft, a mod- 
ification of the basic LQG method was devel- 
oped. 9 *38 A block diagram of the Modified LQG 
design method is shown in figure 6. After a 
state-space model is generated, the first step in 
this approach Is a full -state feedback design. 

Full -state feedback provides for the minimization 
of a quadratic cost function of the output and 
control vectors. To find the optimal full -state 
feedback control law, the quadratic cost function 


Design Methods 


d * E[ X t Q^X + u t Q 2 u] 


( 16 ) 


A flow chart of the overall control law 
design process is shown in figure 5. The first 
element of the process is the selection of design 
objectives (i.e., stability margins, maximum 
allowable structural and control loads, etc.). 

The second element is the selection of a design 
point (i.e., Mach number and altitude). The third 
element is the synthesis of the control law. 

There are several approaches to this key element 
in control law design. Among the approaches are 
classical control theory, modern control theory, 
numerical optimization methods, the method of 
constraints 3 ^, least-squares synthesis 3 ^, and 
the aerodynamic energy method.** The first three 
approaches to this key element will subsequently 
be discussed. The next element, analysis, 
provides information on the performance of the 
control law at off -design conditions. If the 
design objectives are not met, then a gain 
scheduler, which may be a function of Mach number 
and/or dynamic pressure, is evaluated. If use of 
a gain scheduler will still not result in meeting 
the design objectives, the control law synthesis 
element is reentered. 


is minimized. This leads to a control law of the 
form 

u = -C co X (17) 

where C co is the full -state feedback gain 

matrix. 30 Direct measurement of all state vari- 
ables of an aeroelastic system is not feasible. 
Therefore, in the next step, a Kalman estimator is 
used to estimate the state variables from avail- 
able measurements. 

The estimator dynamics are given by 

X c * [ A-BuC c0 - LC] Xc + LY (18) 

where L is the Kalman estimator gain matrix. 30 
However, systems designed using a Kalman estimator 
can have poor gain and phase margins and an 
undesirable high bandwidth. To improve the sta- 
bility margins during the estimator design, the 
"input noise" procedure of Doyle and Stein 3 ^ 
can be used. 


Classical methods . Classical design tech- 
n i ques based on root 1 ocus or Bode plots have been 
applied primarily to single-input-single-output 
(SISO) problems or to multi -input-multi -output 
(MIMO) problems where the coupling between loops 
was sufficiently weak to allow the design to be 
treated as a sequence of single-loop problems. In 
such cases, the classical techniques allow direct 
manipulation of loop gain margins, phase margins, 
etc. Recent classical MIMO design technique 
development includes the characteristic loci 
approach33 and design in a transformed domain 
where input/output pairs are approximately 


Equations (17) and (18) (with X replaced by 
its minimum variance estimate Xq) constitute the 
optimal controller which has the same order as the 
plant model used for the synthesis. As stated 
previously, in the case of a flexible airplane 
model that contains a large number of structural 
modes, unsteady aerodynamic lag states, and actu- 
ator states, the high order of the optimal con- 
troller imposes an unnecessary implementation 
burden. Several sources 9 > 25 > 27 »38,40,41 have- 
shown that a reduced-order controller that 
approximates the full -order optimal controller can 
be found and used with little degradation in the 
closed-loop performance. 
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The transformation of the controller to block 
diagonal form can be used to help select the 
states that are to be retained during the con- 
troller reduction process, A modal residual iza- 
tlon technique 9 can be used to reduce the order 
of the controller. 

The reduced-order controller is then analyzed 
with an evaluation model to examine its perform- 
ance. If the design objectives are not met, then 
there are three paths that can be taken. A 
different set of design variables/controller order 
can be selected, the noise intensities can be 
changed, or the weighting matrices can be changed 
and the optimization and analysis steps repeated. 
The selection of which path to take Is based on 
engineering judgment. 

Numerical Optimization Method . A block 
d i agram of the nume r 1 c a f op t i mi z a 1 1 on method is 
shown In figure 7. The objective of controller 
design is to provide a closed-loop system that 
meets performance specifications. It is natural, 
therefore, that the process be initiated by 
stipulation of the design objectives. These 
objectives can be expressed as inequality 
constraints on variables such as singular values 
of the return difference matrix, control position 
and rate limits, loads, etc. 27,42,43 

The powerful LQG methodology can be employed 
to obtain a candidate stabilizing full-order 
controller. In this phase of the design effort a 
state space model of the plant is employed. It 
should be noted that the plant may contain 
augmented states associated with control logic 
introduced to satisfy part of the design 
criteria. At this point one may choose either to 
refine the design within the MLQG block described 
in the previous section or to exit this block with 
a candidate reduced-order controller. 

The next step is to select parameters in the 
controller which may be employed as design 
variables (e.g., poles and zeroes or coefficients 
of polynomials). One may wish to express the 
controller as a matrix of transfer functions since 
a state space representation is not essential in 
the remaining blocks of this synthesis method. 

Note the alternate path to this point which admits 
the possibility of beginning with a given 
controller form to be specified by the designer. 

The alternate path might be selected to modify an 
existing controller when the plant has changed or 
to develop a scheduling for selected controller 
variables to improve performance at off-design 
conditions. 

The analysis of controller performance for a 
given set of values of the design variables can be 
performed by using both classical and state space 
methods. In the former case, one can avoid the 
s-plane approximation, directly utilize 
experimental frequency response representations of 
dynamic elements, and more readily include the 
evaluation model in robustness assessments. In 
the latter case, one can more reliably obtain 
eigenvalues if they are needed to evaluate 
explicit constraints: the state-space formulation 

may allow more efficient computation of statistical 
or deterministic dynamic response 
characteristics. State-space formulations are 
also more suitable for analytic computations of 


the sensitivities of constraint variables, to 
variations of controller design parameters.. 

The search for a satisfactory controller can 
be automated by using a nonlinear programing 
algorithm 44 * 45 * 46 to determine how to Increment 
the design variables. The search proceeds itera- 
tively and, if a design which satisfies the cri- 
teria is found, the process is terminated. It 
admits explicit consideration of the criteria 
within the design algorithm. If the iterative 
search fails to achieve the design objectives, one 
must try alternate design variables, or modify the 
controller form or relax the design criteria or 
some combination thereof. 

Tool s 

Significant effort at the NASA has been 
applied toward developing tools to help the active 
controls designer in the synthesis and analysis of 
complex active control systems. These efforts 
have resulted in the development, either by the 
NASA or under the NASA sponsorship, of a number of 
computer programs. Three of these programs will 
be described In some depth and several others will 
be identified. 

Analysis Tools 

~ iSAc Tnteractl on of Structures , Aerodynamics 
and Controls . To facilitate the ana lysis of 
active control systems. It is necessary to 
describe numerically, the interaction between 
structural, aerodynamic, and control forces. This 
interaction is described in terms of stability and 
response characteristics. This capability, in the 
form of a computer program system, has been 
assembled and packaged and is identified by the 
acronym ISAC. 47 The system Is in reality an 
assembly of several programs tied together through 
a common data complex as shown in figure 8. 

An analysis proceeds from modal characteris- 
tics obtained separately by any suitable vibration 
analysis program to the aero/structure interface 
(DUN in figure 8) where modal deflections and 
slopes are calculated at points required by a sub- 
sonic doublet lattice code. 48 Generalized aero- 
dynamic forces for purely oscillatory motion are 
then computed as a function of reduced frequency 
and Mach number by using the doublet lattice code 
(DLAT in figure 8). Alternatl vely, both the 
aero/structure interface and the doublet lattice 
code may be bypassed and unsteady aerodynamic 
forces may be input (through the Data Complex 
Manager (DCM) in figure 8) from another source. 

The equations of motion are represented in 
either the frequency domain or state-space 
formulation. Sensor, actuator, and controller 
dynamics can be characterized either in terms of 
transfer matrices or a corresponding state-space 
representation. Four basic types of dynamic 
response (DYNARES in figure 8) analyses can be 
performed within ISAC - stability, stochastic 
response, deterministic response, or frequency 
response. 

An automated capability to determine 
eigenvalue loci as a function of altitude, 
density, velocity, or gain is Included In the 
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stability analysis* If the frequency domain 
formulation Is being employed* the loci are 
obtained either by a determinant or a matrix 
Iteration process with Q-Q(ik) (p-k method) or 
Q*Q(p) (P-P method)* This allows specific roots 
of Interest to be traced* In the state-space 
formulation, a standard eigenvalue problem is 
solved* Zeroes may be calculated by standard root 
determination of the output transfer functions in 
the frequency domain formulation or by solving a 
generalized eigenvalue problem in the state space 
formulation. 

Examples of the stochastic responses obtain- 
able are output rms values for loads, control 
rates and deflections, and sensors due to a unit 
rms input gust velocity or control deflection* 
Power spectral density plots of the outputs also 
may be generated. Output sensitivities to changes 
in specified parameters such as control logic fil- 
ters may be obtained by approximate finite differ- 
encing In the frequency domain formulation* The 
sensitivities may be obtained by closed form 
partial differentiation of steady-state covariance 
performance Indices In the state-space formula- 
tlon.^O 

Output and state time history response due to 
a deterministic input such as a discrete gust or a 
step control input may be obtained. In the fre- 
quency domain formulation, this is accomplished by 
computing the inverse fast Fourier transform of 
the output Fourier transform of interest. In the 
state-space formulation, this can be accomplished 
using a block diagonal form of the transition 
matrix and a convolution integral. 

Performance, stability, and robustness cha- 
racteristics may be determined for both SISO and 
MIMO systems by frequency domain computations. 

Gain and phase margins as well as stability can be 
obtained for the SISO control system from tradi- 
tional Nyquist and/or Bode plots. Analogous MIMO 
stability, performance, and robustness information 
may be determined from certain singular value and 
generalized Nyquist plots. Computations of eigen- 
values and minimum and maximum singular values of 
the return difference and loop gain transfer 
matrices may be obtained from either the frequency 
domain or state-space formulation of the system. 
Frequency responses of these variables as well as 
other outputs of interest may be displayed 
graphically. 

DYLQFLEX - Dynamic Loads of Flexible 
Ai rpl anes wl th Acti ve Control sT DYLOFLEX 49 is a 
system of computer programs which performs 
dynamic loads analyses of flexible airplanes with 
active controls, DYLOFLEX employs the frequency 
domain formulation of the equations of motion. 

The doublet lattice method is employed for 
calculation of the aerodynamic forces. The loads 
equations in DYLOFLEX are developed using the 
method of summation of forces. DYLOFLEX 
incorporates a range of analysis capabilities 
which include calculating dynamic loads due to 
continuous atmospheric turbulence, discrete gusts, 
and discrete control inputs. The output of 
DYLOFLEX consists of statistical quantities of the 
dynamic loads and time histories of the dynamic 
loads. The active control system equations are 


Incorporated by adding rows and columns to the 
basic second-order equations of motion. 

Therefore, all of the active control equations 
have to be expressed as second-order (or less) 
equations, 

NASTRAN - NASA Structural Analysis Program. 
The capability to perform analyses of an actively- 
controlled flexible airplane was added to the 
NASTRAN program in 1979. 29 NASTRAN employs the 
frequency domain formulation of the equations of 
motion. Four different methods are available for 
the calculation of the aerodynamic forces* These 
include the doublet lattice method, strip theory 
method, Mach box method, and piston theory 
method. The different analyses available include 
flutter analyses, transient response, and random 
response analyses. Active control system 
equations are incorporated in a similar manner as 
DYLOFLEX where they must be expressed as 
second-order (or less) equations. 

Synthesis Tools 

SYNpAC/PADLQCS - Synthesis Programs for the 
Design of Active Controls 

SYNPAC - Synthesis Package for Active 
Controls 

PADLOCS - Program for the Analysis and Design 
of Linear Optimal Control Systems 

PADLOCS and SYNPAC are two packages of 
programs which have been developed at the 
NASA-Langley to facilitate the synthesis of 
control laws. Both the MLQG and constrained 
optimization techniques described previously are 
available in these programs. Both PADLOCS and 
SYNPAC rely heavily on ORACLS 50 software for 
the LQG portion of the design. Both systems of 
programs can be used to perform nonlinear 
constrained optimization studies to seek an 
implementable reduced-order controller which meets 
a set of design criteria. The design criteria are 
included as inequality constraints. Various 
parameters such as rms control rates and minimum 
singular values may be included in a weighted 
performance index. 

The SYNPAC system makes direct use of I SAC as 
shown in figure 9. SYNPAC and I SAC share a 
common data base and data base manager, and the 
DYNARES module of I SAC is employed to perform the 
analyses. Consequently, the constrained 
optimization portion of SYNPAC can utilize either 
a frequency domain or state-space representation 
of the plant and controller. Stability can be 
constrained within SYNPAC either by eigenvalue 
computations or by computing the number of 
encirclements of the -1 point in a generalized 
Nyquist diagram. Also in SYNPAC, constraints on 
parameters such as short period damping may be 
included. The performance index Is a weighted sum 
of the squares of the magnitudes of selected 
outputs, augmented by one of several types of 
penalty functions which incorporate the constraint 
violations. The user may select from a non- 
gradient 51 or several conjugate gradient based 
nonlinear programing algorithms 45 to perform the 
constrained optimization. In actual application, 
only the nongradient based algorithm has been 
employed. Consequently , no effort has been 
expended to develop efficient code for gradient 
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computations; currently, gradients are computed 
numerically by finite differences. SYNPAC has 
been employed in several control law design 
studies. 25 ** 7 


A flow chart for the PADLOCS program is shown 
in figure 10.. PADLOCS requires a state-space 
representation of the plant and controller. Plant 
matrices ([A],[B U , B wg U and [C] in eq. 8), 
which could be developed by ISAC are accepted as 
input to PADLOCS. A feasible directions, conju- 
gate gradient algorithm is employed by PADLOCS to 
optimize a candidate reduced-order controller. 
Closed-form expressions for the gradient informa- 
tion are coded in PADLOCS for constraint and func- 


tion variables involving output rms and return 
difference singular values. Stability is examined 
at each iterative step by solving for the 
eigenvalues. The PADLOCS system has been employed 
to design a number of active control laws. 


Comparison of Analysis With Experiment 

Active control experiments (wind tunnel and 
flight) have been used to validate analysis, to 
evaluate feasibility, and to demonstrate predicted 
benefits. A description of the wind tunnel and 
flight experiences to date with active controls is 
given in reference 52. The majority of the 
wind-tunnel experiments have focused on flutter 
suppression. Flight experiments have concentrated 
on the other active control concepts. A few 
representative examples of these tests will be 
described where analytical results have been com- 
pared with experimental data. 


control system on the C-5A fleet designed to 
reduce wing loads. 

A description of both the analytical and 
flight test data for the DAST vechicle is given in 
reference 54. A comparison of both frequency and 
damping of the dominant mode is shown In figure 
13. The analysis and experimental data are for an 
altitude of 15,000 feet. The change in frequency 
with Mach number is predicted well for both the 
FSS-off and FSS-on cases. However, analysis 
overpredicts the damping for both the FSS-off and 
FSS-on cases. The experimental flutter speed is 
extrapolated to be approximately M » 0.80 for the 
FSS-off case. An actual flutter point was 
encountered for the FSS-on case at H « 0.82. The 
analysis overpredicts the FSS-off flutter speed by 
4 percent and overpredicts the FSS-on flutter 
speed by 2 percent. 

A description of the C-5A ALDCS is given In 
reference 15. Two examples of the comparison 
between analysis and experimental data are taken 
from reference 15. figure 14a shows the amplitude 
of wing root bending moment frequency response 
function. The flight condition is 300 knots at an 
altitude of 25,000 feet. There is good agreement 
between analysis and experiment for both ALDCS-off 
and ALDCS-on cases. Figure 14b shows a comparison 
between analysis and experimental data for the 
spanwise distribution of incremental maneuver 
bending moment. The flight condition is a Mach 
number of 0.78 and an altitude of 30,000 feet. 
Again, there is good agreement between the anal- 
ysis and experimental data. 


Wind tunnel . 

The - wind-tunnel example is a DC-10 derivative 
model that is described in reference 53, The 
model was designed to be tested in a low speed 
wind tunnel and therefore no transonic effects 
were present. One of the purposes of the tests 
was to assess the accuracy of dynamic analysis 
methods applied to the active control functions of 
flutter suppression and gust load alleviation. 
During the tests, parametric variations of control 
law gain and phase were made and experimental data 
were gathered. Figure 11 shows a comparison 
between the experimental and predicted stability 
boundaries as a function of gain and phase. 

Figure 12 shows a comparison between experimental 
and predicted gust loads as a function of gain and 
phase. The stability boundary predictions 
employed the state-space formulation for the equa- 
tions of motion whereas the gust loads predictions 
employed the frequency domain formulation. The 
comparison of predicted and experimental results 
indicated that all stability and gust loads trends 
(in terms of gain, phase, and velocity) were 
predicted properly and that the analysis gave 
realistic estimates of levels. 

Flight 

‘" Two flight tests of active control concepts 
that involve significant structural response will 
be discussed. The first is NASA’s DAST (Drones 
for Aerodynamic and Structural Testing) program 
which is a research program whose primary 
objective is to develop and evaluate active 
control analysis and synthesis techniques. The 
second is the ALDCS (Active Load Distribution 
Control System) which is now an operational active 


Future Trends 

The following is a brief discussion of the 
authors' perspective on future trends In analyti- 
cal methods. These items are by no means exhaus- 
tive but are only considered to be 
representative. 

Since the critical speed range for flutter is 
often transonic and the use of supercritical air- 
foils is increasing, the Interfacing of the new 
nonlinear transonic unsteady aerodynamic force 
computations with structural dynamics and control 
system equations will become increasingly 
important. 55 The interfaces will be needed for 
both analysis and control law synthesis. In the 
control system design area, future emphasis will 
be on functionally integrated control system 
design for total mission control. In the robust 
control law synthesis area, the generalization of 
the structured singular value problem 55 , optimal 
projection 57 and constrained optimization meth- 
ods for fixed-order control law synthesis for 
both continuous and sampled-data systems will be 
gaining more attention. 

Concluding Remarks 

Analytical methods employed in active 
controls analyses and design have been described. 
Development of the math models from the basic 
frequency domain equations of motion to the state- 
space equations of motion have been discussed. 

Both classical and modern control techniques, as 
applied to the active controls problem, have been 
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reviewed, A combination of these techniques has 
been receiving considerable attention recently. 
Tools that Implement these techniques have been 
described. Future trends toward functionally 
integrated control design, sampled-data 
techniques, robust control law synthesis 
techniques, and the incorporation of transonic 
unsteady aerodynamics have been identified. 
Examples where analytical results have been com- 
pared with experimental data have been presented. 
For these examples, the comparisons have been 
quite good. However, for these examples, the 
math models employed for the control law design 
have had the benefit of previous experimental 
results. The advances that have ocurred in analy- 
tical methods are allowing certain nonflight -cri- 
tical active control functions to be considered 
during aircraft preliminary design. However, con- 
siderable work remains to be done before math 
models can be used with the necessary confidence 
for flight-critical active control functions. 
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FIGURE 4. Comparison of p-k and p-p Techniques 
for Determining Stability 



FIGURE 5, Flow Chart of Overall Control Law 
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FIGURE 8. Flow Chart of ISAC Analysis Tool 
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